We define a block persistence probability pi(t) as the probability that the order parameter 
integrated on a block of linear size I has never changed sign since the initial time in a phase ordering 
process at finite temperature T < T c . We argue that pi(t) ~ l~ zB ° f(t/l z ) in the scaling limit of 
large blocks, where 80 is the global (magnetization) persistence exponent and f(x) decays with the 
local (single spin) exponent 9 for large x. This scaling is demonstrated at zero temperature for the 
diffusion equation and the large n model, and generically it can be used to determine easily 80 from 
simulations of coarsening models. We also argue that 80 and the scaling function do not depend 
on temperature, leading to a definition of 8 at finite temperature, whereas the local persistence 
probability decays exponentially due to thermal fluctuations. We also discuss conserved models for 
which different scaling are shown to arise depending on the value of the autocorrelation exponent 
A. We illustrate our discussion by extensive numerical results. We also comment on the relation 
between this method and an alternative definition of 8 at finite temperature recently introduced by 
Derrida [Phys. Rev. E 55, 3705 (1997)]. 
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I. INTRODUCTION 

Phase ordering processes correspond to the dynam- 
ics of systems quenched from a disordered high temper- 
ature state to a temperature where the equilibrium state 
is ordered. Dynamics proceed through coarsening of or- 
dered domains, and the domain linear scale L(t) diverges 
as t 1 / 2 . In the coarsening regime, nontrivial spatial and 
temporal correlations develop, adopting a scaling form. 
For instance, the order parameter equal time correlation 
function (^(x,iM(M)) = /(|x|/L(t)). 

Universality classes depend not only on the space di- 
mension and the symmetries of the order parameter, as 
for static critical phenomena, but also on the conserva- 
tion laws of the dynamics. Indeed, for a single static 
universality class, several dynamics can be used with 
the only constraint that they must obey detailed bal- 
ance. For a scalar order parameter, non conserved dy- 
namics (model A) describe a ferromagnet, while con- 
served dynamics (model B) describe demixtion or seg- 
regation in binary alloys. Consequently, the set of dy- 
namical critical exponents such as z and A, defined by 
(</j(x, t')<p(x, t)) ~ [L{t) / L{t')] x , for t' > t, is not related 
to static exponents by any hyperscaling law. 

A remarkable point is that as far as the temperature T 
of the quench is concerned, there are only two universal- 
ity classes, namely T = T C (critical quenches) or T < T c . 
This is assessed by numerics, renormalization group re- 
sults or large-n expansions Q . More precisely, two point 
correlations have the same scaling (up to multiplicative 
constants) for any T < T c . Therefore the temperature 
is an irrelevant parameter for T < T c quenches, and one 
may set T = as well. 

However, the situation is not as simple if one considers 
quantities involving more subtle correlations. One such 
quantity, which has attracted much interest recently, is 
the persistence probability, defined as the probability that 
the local order parameter at a given point x has never 
changed sign since the initial time ^-^]. For instance, in 
simulations of the Glauber Ising model, p(t) is the frac- 
tion of spins that have never flipped since the initial time. 
At T = 0, p(t) is the probability that a given point has 
never been crossed by a domain wall. It usually decays 
with an exponent p(t) cx t~ e . For general nonequilib- 
rium dynamics, p(t) is the probability that a zero-mean 
stochastic quantity has never changed sign since the ini- 
tial time. The analytical study of p(t) is difficult due 
to the fact that it probes the whole history of the pro- 
cess. Even for simple scalar diffusion with a zero-mean 
random initial condition, a nontrivial algebraic decay is 



found At T c , the persistence of the global magne- 

tization was shown to yield a new independent critical 
exponent for the Ising model j^U ■ 

The temperature universality of the T < T c scaling of 
correlations, corresponding to a single fixed point T = 
in the renormalization group, seems to be broken for the 
local persistence p(t), because at T > 0, thermal fluctua- 
tions lead to an exponential decay of p, in contrast to the 
power law decay at T = 0. To address this question, Der- 
rida || recently proposed to study persistence at finite 
temperature for nonconserved Ising and Potts models by 
comparing two systems A and B evolving with the same 
thermal noise from two different initial conditions. Sys- 
tem A is initially in a completely random configuration 
whereas B is in its fundamental (all spins assuming the 
same value). Persistence is now defined as the proba- 
bility r(t) that SfSf(t) has kept a constant sign since 
t = 0. The underlying idea is to discard simultaneous 
flips (at the same site) in both systems, because flips in 
B are only due to thermal fluctuations. 

The implementation of this simple idea by Derrida || , 
and more extensive simulations performed by Stauffcr 
p0| , have shown that r(t) decays algebraically. The ob- 
served exponent seems to be temperature independent 
and equal to the T — local persistence exponent 9, for 
the Ising model. Therefore, universality seems to hold 
with this new definition of persistence. However, the 
method cannot be used for conserved models, as system 
B would not evolve. Since Kawasaki (spin-exchange) dy- 
namics freeze at zero temperature, a definition of persis- 
tence at finite temperature applying to conserved models 
is required. Moreover, Derrida's definition is not easy to 
generalize to continuous models. A definition involving 
a single system would be more satisfactory, as we know 
from the study of damage spreading that behaviors of ob- 
servables obtained by comparison of two systems evolv- 
ing with the same noise often depend on the Monte-Carlo 
algorithm used (see below). 

In a recent letter [jll), we introduced the notion of 
block persistence as a very natural method to give a 
temperature independent and intrinsic definition of the 
persistence exponent. The method is in a way an d la 
Kadanoff implementation of the renormalization group 
ideas underlying the universality of correlations. The 
block persistence probability pi (t) is the standard persis- 
tence probability for a coarse-grained variable obtained 
by integrating the order parameter on a block of linear 
size I. In fPH , we argued that the large / scaling of pi(t) 
is independent of T and corresponds to the T = fixed 
point, because increasing I reduces the relative thermal 
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fluctuations of the block variables. 

In this article, we give a more detailed and general 
discussion of block persistence, which we illustrate with 
extensive simulations of different coarsening models. The 
structure of the paper is the following. We start by re- 
viewing in section O a few mathematical results needed 
to discuss persistence for physical models. 

In sec. Ill, we comment further on Dcrrida's compar- 
ison method and check its intrinsicality. We show that 
even if the persistence exponent does not seem to de- 
pend on the algorithm used, the cross-over to T c and the 
T > T c behavior of the persistence probability is com- 
pletely different for heat bath and Glauber dynamics. 



In sec. IV, we start from a general discussion of global 



persistence below T c , and define block persistence as a 
natural way to include in a single framework the global 
and the local persistence exponent, through its scaling for 
I oo with l/L(t) fixed at T = and T > 0. At T = 0, 
we explicitly prove the postulated scaling form for the 
diffusion equation and the large-n model. We show that 
block scaling leads to an easy numerical determination of 
the global persistence exponent 9q. We present numer- 
ical results for several systems, illustrating the previous 
discussion. 

In Sec. we move to finite temperature and justify 
that the scaling should be the same as at zero tempera- 
ture, because the thermal exponential decay is eliminated 
in the scaling limit of large blocks. Thus block persis- 
tence provides with a definition of local persistence at 
finite temperature. We present simulations for the Ising 
and Potts models, illustrating temperature universality. 
We also discuss the T = T C case. 

In Sec. VI, we discuss the special case of conserved 
order parameter dynamics. Block scaling works as for 
the nonconserved case, but for an important feature: we 
analytically predict that the scaling function should be 
qualitatively different for A = d and A < d systems. This 
prediction is confirmed by simulations of one-dimensional 
models. We also present finite temperature simulations 
for the two-dimensional Kawasaki dynamics. 



II. MATHEMATICAL AND GENERAL RESULTS 

Before moving to physical problems, we would like to 
summarize a few useful mathematical results. Consider 
a general stochastic process X(t), with (X(t)) = 0. We 
are interested in the probability p(t) that X(t') > for 
all < t < t'. 

This is an old problem in probability theory fll2| , |l3"| | , 
but a difficult one, and despite the large number of papers 
devoted to this subject, very few quantitative results are 
known, most of them concerning stationary and Gaus- 
sian processes, which are completely determined by their 
correlator C(r) = (X(t)X(t + r)). With these strong 
restrictions, p(t) still cannot be computed analytically, 
even in the large t limit. Actually, p[t, C(r)] is known 



only for very few specific correlators [ |12[|13| ] . One of these 
correlators is C(r) = e~ ar , which is the general correla- 
tor of a Markovian stationary Gaussian process with the 
condition C(0) = I, for which 



p(t) = - arcsin(e at ) 

7T 



(1) 



and p{t) ~ (2/7r)e" at at large t. 

Generally speaking, p(t) and its asymptotic large t de- 
cay depend sensitively on the whole function C(r) and 
not only on its behavior for small or large r. For instance, 
Majumdar and Sire have considered the Gaussian 
process with C(r) = (1 — e)e~ T + ee~ 2r . Despite the fact 
that C decays oc e~ T at large r for all e < 1, simula- 
tion of the process shows that p{t) oc exp(— a(e)t), where 
a(e) interpolates continuously from 1 to 2 when e is var- 
ied from to I. In (li), the Markovian correlator was 
used as a starting point for perturbative and variational 
approximations, which are however uncontrolled. 

The following rigorous results for any stationary Gaus- 
sian process with zero mean are also very useful [ 12l 1J 



p[t,bC(r)} = p[t,C(r)} (2) 
p[t,C(bT)} = p[bt,C(r)} (3) 
(Vr, d(r) > C 2 (r)) (Vt, p[t, Ci(t)] > p[t, C 2 (r)}) (4) 

From the first relation, we see that p is completely de- 
termined by the normalized correlator C(r)/C(0). The 
second relation will be used to obtain scaling forms for 
persistence probabilities in the following. Finally, the 
third relation shows that p(t) decays exponentially in 
time for a stationary Gaussian process with a correla- 
tor that is bracketed for all r by two Markovian corre- 
lators e~ b ' r ' < C(t) < e~ a l r l, because then p(t) is also 
bracketed by two exponentials. Most of the correlators 
encountered in physical noncquilibrium processes actu- 
ally have this property in a proper time variable (see 
below). However, there might be power law prefactors in 
the large t decay of p(t) . 

Of course, in nonequilibrium dynamics, stochastic pro- 
cesses are scarcely Gaussian, and, by definition, never 
stationary in physical time. However, if there is scaling 
relatively to a dynamically diverging scale L(t), one must 
have for large t and t' 



a(t,t') 



(X(t))X(t>)) 



f[L(t)/L(t% (5) 



with f(x) = f(l/x). This implies the stationarity of the 
process X(t)/ \J (X 2 (t)) in the variable u = \aL{t). 

Now if the process is Gaussian, we obtain that gener- 
ically p(u) decays as e~ 9u and therefore p(t) decays as 
L(t)~ e . For most systems L(t) oc t x ' z and we recover the 
power law decay in time with 9 = z9. The simplest exam- 
ple of such a Gaussian process is the diffusion equation 
(see below). Still, because the process is non Markovian, 
9 cannot be computed analytically, and an independent 
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interval approximation was used to predict accurately 9 

To end with this general discussion, we consider the fol- 
lowing situation, which will be of use in the study of block 
persistence. Consider a family of Gaussian processes in- 
dexed by a variable I > 0, {Xi(t)}, with normalized cor- 
relators ai(t,t'), with the following scaling property 



Sf (t + Ai)=sign 



l + tanhO^SfW) 



(8) 



ai(t,t') = h{t/l z ,t'/l z 



(6) 



Then obviously Xi(t) 
Pi(t/l z ). 



X%(t/l z ), leading to p(t) 



III. COMPARISON OF SYSTEMS 

Now, let us come back to coarsening processes. Con- 
sider the nonconserved Ising dynamics. At T = 0, 
p(t) cx t~ s , where 8 is nontrivial and seems to be in- 
dependent of other exponents. This is due to the fact 
that spins cannot flip when they are within an ordered 
domain. Flips occur only at interfaces between domains, 
and the slow surface tension driven motion of these in- 
terfaces makes for the slow decay of p. 

The situation is dramatically different at finite temper- 
ature, because thermal fluctuations allow energetically 
forbidden flips. These activated flips occur with a de- 
cay rate r ~ e -AE/k B T ^ ^gj-g is a typical energy 
barrier to flip a spin inside a domain, of order the ex- 
change constant J. Therefore, these thermal flips lead to 
an exponentially decaying p(t) cx e~*/ r (see sec. 0). 

Why then is there a unique scaling of correlations at 
finite T < T c ? The reason is that the domain struc- 
ture in the scaling regime is the same at any T < T c . 
The thermal fluctuations cancel out in the two point cor- 
relations, which reflect only the alternation of domains 
of different phases. The temperature dependence of the 
value of the bulk magnetization (approximately equal to 
its equilibrium value), just leads to a temperature depen- 
dent multiplicative constant in the scaling function. 

From this point, it becomes clear that a simple tem- 
perature independent definition of 9 should be through 
the probability r(t) that a given site has never changed 
phase, i.e. has never been crossed by a domain wall. At 
T = 0, we clearly have r(t) = p(t), and at T < T c because 
of the universality of the domain dynamics, r(t) should 
have the same decay as at T = 0. 

Derrida J^] proposed a very clever scheme to implement 
this idea for the nonconserved Ising model by simulating 
two systems A and B evolving with the same Monte- 
Carlo dynamics, with the same thermal noise. System 
A is prepared in a completely random initial condition, 
whereas B is prepared in the fundamental state (all spins 
equal to one). Then both systems are updated simul- 
taneously using the heat bath algorithm with the same 
random number z at the same site i: 



Sf(t + At) = sign 



l + tatjhG9£i W)) 



(7) 



Then, the fraction of persistent spins r(t) is defined 
as the fraction of sites for which SfSf has kept a con- 
stant sign since t = 0. It means that we discard flips 
that occur simultaneously in both systems, because flips 
in system B are purely thermal fluctuations, as there is a 
single -I- phase. Accordingly, Derrida found that at finite 
temperature T < T C) r(t) cx t~ e , with 8 consistent with 
the T — persistence exponent. This was confirmed by 
extensive simulations performed by Stauffer jnj . 

However, this practical definition of persistence for the 
Ising model is not completely satisfactory. First, it can- 
not be directly adapted to continuous models. Indeed, 
for a continuous order parameter the probability of a si- 
multaneous flip in both systems will be zero in continuous 
time. 

A further restriction is that the method cannot be used 
for conserved dynamics, as the Kawasaki spin-exchange 
dynamics, because system B would not evolve from a uni- 
form initial condition. There is no proper initial condi- 
tion for system B. This restriction is important, because 
Kawasaki dynamics cannot be studied at zero temper- 
ature, and therefore a definition of persistence at finite 
temperature is required. 

Finally, one would be more satisfied to get an intrinsic 
definition of persistence. The comparison of two systems 
evolving from different initial conditions has attracted 
much attention, especially in relation to the notion of 
damage spreading . It was soon realized that the be- 
havior observed depends on the implementation of the 
Monte-Carlo algorithm. Therefore, one could fear that 
Derrida's definition may work only with the heat bath 
algorithm. To check this, we performed simulations with 
the heat bath algorithm and the Glauber algorithm. We 
find that the T < T c behavior is the same for both dy- 
namics. However, quite interestingly, the T > T c behav- 
ior of r(t) is completely different. 

For heat bath dynamics, for T > T c , r(t) reaches a 
plateau. This corresponds to the fact observed by Der- 
rida and Weisbuch that above T c two systems evolv- 
ing with this algorithm become identical within a finite 
time. When T — > T c , this plateau crosses over to a power 
law r(t) ~ t~ 6c . From simulations at T c we find 9 C « 0.9, 
but 9 C can also be extracted from a scaling analysis of 
the cross-over for T -> T+ . At finite T > T c , there 
are no domain walls. Starting from an infinite temper- 
ature state with a correlation length £ = 0, £ increases 
to reach its equilibrium value £ eg . In the vicinity of T c , 
~ (T — T c )~" is very large. Therefore, at early times, 
for £(t) <C (T — T c )~ v the system behaves as if it were to 
reach a critical (infinite £) equilibrium state, i.e. as if it 
were at T c , and £(t) ~ t 1 ^" while r(t) ~ t~ 6c . Deviations 
from this power law behavior appear only at late times 
when £(t) approaches the finite value £ e9 and r(t) reaches 
a plateau. Consequently we expect the scaling form 
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r(t) ~ {^ q ) a gUt)IU ~ t" e °/[t(T - Tc)«" 



(9) 



where /(a;) cx x e<: when x — + oo and f{x) tends to a 
constant when x — > 0. 
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FIG. 1. Scaling behavior of the fraction of persistence spins 
r(t), from Derrida's definition, for the Ising model with heat 
bath dynamics when T — » T c . Simulations were carried out 
on a 1000 2 lattice and 20 samples were averaged. The data 
collapse is obtained with 8 C — 0.9, a = vz c , with the exact 
value v = 1, and z c — 2.17. The scaling function goes to a 
constant at small argument and diverges as a power law at 
large argument. 



10" 



10" 



■ :o- 



1(T 




10 



100 



1000 



FIG. 2. Decay of the fraction of persistent spins r(t), with 
Derrida's definition, for the Ising model with Glauber dynam- 
ics at T — T e , 1.005T C and 1.1T C , from simulations on a 1500 2 
lattice (20 samples). Above T c , r(t) decays faster than any 
power law, in contrast with the heat bath dynamics. At T c , 
r(t) also seems to decay faster than algebraically, but one 
cannot positively rule out a power law decay with a large ex- 
ponent 6 C bigger than 2.4 (the extrapolated exponent from 
our data). 



This scaling behavior is illustrated in figure EL which 
shows results of simulations of heat bath dynamics for the 
two dimensional Ising model at different temperatures 
above T c . The best scaling is obtained with 9 C = 0.9 
(v = 1 and z c = 2.17). It is quite surprising to ob- 
tain a new exponent at T c (see the discussion for block 
persistence below), and one should wonder whether this 
exponent is universal or if it depends on the chosen algo- 
rithm. 

If we consider another frequently used algorithm, the 
Glauber dynamics which corresponds to the update rule, 

"1 + taring £\ ,%(*)) 



S t (t)=Si(t) x sign 



(10) 



and using the same z at the same site for system A and 
B, we find a very different behavior of r(t) above T c . As 
shown on figure |^, it now decays faster than any power 
law. At T c , it is difficult to distinguish from our simu- 
lations performed on a 1500 2 lattice, averaging over 20 
samples, whether r(t) decays exponentially or as a power 
law with an exponent 9 C bigger than 2.4 (the value ex- 
trapolated from our data). Anyway, we do not find the 
exponent 9 C = 0.9 found for heat bath dynamics. Thus, 
this exponent is not intrinsic, neither is the T > T c be- 
havior of r(t) (similar results have also been found by 
Hinrichsen and Antoni Jl7|). 

This illustrates the kind of problems that can be en- 
countered using observables defined by the comparison of 
two systems. Note however that, as said before, the large 
t decay of r(t) below T c is the same for both dynamics, 
and therefore seems to be intrinsic. The cross-over in the 
vicinity of T~ will be, of course, different. One could be 
tempted to relate the different behavior obtained above 
T c to damage spreading properties of the dynamics. The 
question of damage spreading is to know whether the 
distance (in configuration space) of two copies of a same 
system, evolving from two slightly different initial con- 
ditions, diverges (damage is said to spread), or keeps 
bounded (damage is said to heal). For the 2d Ising model, 
damage heals for heat bath dynamics and spreads for 
Glauber dynamics. However, in Id, damage heals for 
both dynamics, and we have checked that both dynam- 
ics lead to an exponential decay of r(t). Therefore, the 
absence of damage spreading is not a sufficient condition 
to obtain the saturation of r(t), which seems to be due 
to a very specific property of heat bath dynamics in 2d. 

In the rest of the article, we describe a completely 
different approach to finite temperature persistence, in- 
volving a single system, which therefore avoids such dif- 
ficulties and can be applied to continuous and conserved 
models. 



IV. BLOCK PERSISTENCE AT ZERO 
TEMPERATURE 

We consider the nonequilibrium dynamics of a non- 
conserved order parameter ip(x.,t), which can be cither 
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continuous or discrete. The case of a locally conserved 
order parameter will be postponed to section VI. We 
first discuss block persistence at T = 0, because as we 
will see later the scaling properties of pi are temperature 
independent. 

Actually, introducing block persistence is also quite 
natural at zero temperature, because it provides a link 
between the global and the local persistence probabil- 
ities. Block scaling will be shown below to be a very 
effective way of determining the global persistence ex- 
ponent #0j corresponding to the global order parameter 
M(t) = J (p(x,t)cbc. This global persistence exponent 
has been extensively studied at T c J7|p p^ , pj| . Here, we 
would like to thoroughly discuss the T < T c case. We 
shall assume T = 0, but the discussion would be the 
same at any T <T C . 



A. Global persistence 

First, we remark that while the stochastic process 
{y(x, t)} (at a given point) is generally speaking both non 
Gaussian and non Markovian, {M(t)} is always Gaus- 
sian in the thermodynamic limit. Indeed, the vector 
(M(ti), .., M{t n )) is the sum of an infinite number of vari- 
ables (y(x, ii),..., (p(x,t n )). Since the correlation length 
for ip is finite at t\, ■-,t n , the ip vectors have short range 
correlations and the central limit theorem entails that 
the magnetization vector is Gaussian, for every choice 
of an arbitrary number n of times. This proves that the 
whole process {M(t)} is Gaussian (which is stronger than 
just saying that M(t) is Gaussian at any t). For a finite 
system, there are non Gaussian corrections due to the 
fact that the number of independent contributions to the 
magnetization vector is finite and of order V/L(t n ) d , if 
t n is the largest of the t{. 

Therefore, for an infinite system, the global persistence 
probability is completely determined by the two-times 
correlator 



a g (t,t>) = (M(t)M(0)/#Wf)) (11) 

The analytical determination of 9q is consequently sim- 
pler in principle than for 9 (for 9, a nonlinear Gaussian 
approximation was used by Majumdar and Sire H]). At 
T = T c , Majumdar et al |Q have been able to compute 
an e expansion of the the global exponent 9 C for model A. 
For T < T c , there is no natural perturbation parameter 
as e. The analytical study of 9q can be performed using 
the methods of Q in dimensions d > 3 [ p0| . 

Interestingly there is a relation between the autocor- 
relation exponent A and 9o when M is a Markov process 



%z = A - d/2. 



(12) 



This relation is the consequence of the scaling of corre- 
lations and is the counterpart of a similar scaling law at 



To show Eq. (|l2|), we use the fact that the Gaussian 
process M(t) / M(i) is stationary in the scaling limit as 
a function of u = \nL(t) (see sec. |l|), with 



a g (t,t')=f[L/L'}=c(\u-u , \). 



(13) 



The two point correlator C(k,t,tf) = L x L /d ^ x g(kL') in 
the scaling regime for t! ^S> t. If g(0) = 0(1), which is the 
case for nonconserved models, (M(t)M(t')) oc [L/L'] x L' d 
while (M 2 (t)} oc L d , and we obtain, 



a g (t,t') 



' L ' 


A-d/2 


V 


(\-d/2)/z 






t' _ 





for t' > t (14) 



Up to this point, the results are general and valid for 
a non Markovian process. Now, if the normalized M is 
Markovian, then necessarily c(\u — u'\) — exp(— 9qz\u — 
u'\) (see sec. ||). In other words f(x) — x~ 9 ° z for all 
x > 1, and since Eq. (yjj) expresses that f(x) ~ a;( d / 2_A ) 
for x 3> 1, we obtain Eq. (p^). Note that the lower bound 
A > d/2 proposed by Fisher and Huse ensures that 
0o is nonnegative. Below, we shall demonstrate that M is 
Markovian for the T = one dimensional Glauber model, 
and we shall find that 9$z = 1/2 = A — d/2 (A = 1), but, 
generally speaking, Eq. (|l2| ) is violated because M{t) is 
non Markovian and f(x) is not a pure power law. The 
Markovian value of #o is neither an u pper nor a lower 
bound (see numerical results in section fVE). The only 
general bound we have is 



< 



(15) 



since obviously to flip the magnetization one has to flip 
single spins. 

The direct determination of 9q is quite difficult. One 
has to record the time when the magnetization first 
changes sign, for a very large number of runs, which limits 
drastically the sample size. Cornell and Sire |p2[ simu- 
lated the two dimensional Ising model on a L — 8 to 128 
lattice, and were obliged to use a finite-size scaling anal- 
ysis that did not prove very conclusive, leading to a large 
uncertainty on the value of 9q « 0.06 ~ 0.11. We shall 
see below that block persistence, which we now define, 
leads to a much easier determination of 9q ■ 



B. Block persistence 

The idea is to define a more general quantity, the block 
persistence probability, that coincides with the global and 
the local persistence in different limits. The procedure 
is very natural: we consider a coarse-grained variable 
(fl(x,t), obtained by integrating scales smaller than I. 
The simplest procedure is to integrate tp over a block of 
linear size I, as will be done for numerical simulations 
of lattice spin models. Alternatively, one can also elimi- 
nate Fourier modes of wavelength smaller than I, as will 
be more convenient for the analytical treatment of con- 
tinuous models. The block persistence probability pi(t) 
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is just the persistence probability for the coarse-grained 
variable. For I — oo we recover the global persistence, 
while for I = (or 1 on a lattice), we get the local per- 
sistence. 

Now, for finite I, the time dependence of pi interpo- 
lates between the two exponents 9 and 8q. Indeed, at 
early times, when L(t) <C I, the system effectively sees 
infinite blocks, and pi(t) oc t~ e " . Then for L{t) ^> I, 
blocks behave as single spins, and pi(t) ~ cit~ e . There- 
fore, we expect a scaling form of pi(t) for I — > oo with a 
fixed ratio l/L(t) 



Pl (t)~r a g (L(t)/i) = r a f(t/i z ), 



(16) 



where f(x) oc x~ e ° when x — » and fix) oc x~~ e when 
x — > oo. a must be equal to z9q because for finite t and 
I — * oo, pi(t) must tend to the global persistence proba- 
bility. 

This scaling form can be demonstrated for two analyt- 
ically tractable models closely related, namely the diffu- 
sion equation and the large- n limit of the O(n) noncon- 
served model. The reason is that in both models, the 
process {ip{x,t)} is Gaussian, entailing that all coarse- 
grained variables are also Gaussian, and pi{t) depends 
solely on the normalized correlator 



ai{t,t') 



( W (x,i)^(x,i')) 



(17) 



which can be computed analytically. 

C. Diffusion equation 

The diffusion equation may be the simplest example 
of nonequilibrium dynamics. It is not really a coarsening 
model, because of the absence of domain walls, due to 
the linearity of the equation. Consider a scalar field tp 
evolving according to 



Tt =v > ' 



(18) 



from a random initial condition with zero mean (tp) = 
and short range correlations (tp(x,0)p(x' ,0)) — A5(x — 
x 1 ). For this model, the global magnetization is con- 
served, leading to 9o = 0. 

Integrating the equation in Fourier space, we obtain 
the Fourier transform of the correlator 

C(k,M') = (0(M)c3(-M)) = Ae~ fc2 (*+*'). (19) 

Computing the two time correlator C(t,t') = 
5^ k C(k, i,t') leads to the normalized correlator 



a(t,t') 



C(t,t') 



Att' 



y/C(t,t)C{t',V) \(t + t>) 2 



(20) 



This correlator yields a nontrivial persistence exponent, 
which can be reproduced with an excellent precision us- 
ing an independent interval approximation (IIA) |||| 

As remarked before, considering blocks of size 
I is equivalent to introduce an upper cut-off in 
Fourier space A ~ l/l, and to consider ipi — 
(1/VV) X)|k|<A ^(k) exp(ik.x). The correlator of the 
corresponding block variables is, 

d (t, = (p t (x, t) W (x, f)> = £ <?(k, i, (21) 

|k|<A 

The ipi variables arc Gaussian, and the behavior of 
pi(t) depends only on the normalized correlator ai(t,t'). 



a\(t,t') = 



VtF \ 2 F(\ 2 (t + t')) 

t + t' y/ F(\ 2 t)F(\ 2 t') 



H(\% A t') 



(22) 



with F(x) = J* y d 1 e y2 dy. From the final discussion 
of sec. ||, we have pi(t) — pi{t/l 2 ), which is precisely 
the scaling form assumed from physical arguments, with 
a = 0. 

The probability p\ (t) cannot be explicited, but we can 
obtain its asymptotic behavior. For t,t' ^> I 2 , or for 
A — > oo one has 



a x {t,t') - a(t,t') 



(23) 



i.e. we recover the one point two-time normalized corre- 
lator, leading to a nontrivial exponent 9. Then, in the 
opposite limit of large blocks (or small times) t,t' <C I 2 , 
one has a\(t,t') = l + 0(X 2 (t + t')), which corresponds to 
a non-evolving field. The scaling of the block persistence 
probability is therefore conform to the general discussion 
above. 



D. Large n limit 

Now, let us investigate the 0{n) model in the large n 
limit. As usual, we start from a n-components vectorial 
order parameter tp with the Time Dependent Ginzburg- 
Landau dynamics, 



d t <Pa = V 2 v?o 



rp a 



9 -2 
-tp a tp 

n 



(24) 



In the large n limit, tp 2 /n can be replaced by the average 
(tp 2 ), where tp is now any component of the field, and one 
obtains a linear self-consistent equation, which reads in 
Fourier space, 



d t p(k,t) = -{k 2 + R(t))tp(k,t), 



(25) 



with R(t) = r + g(tp 2 ). 

Hence, <p(k,t) = tp(k, 0)/i(t)" 1/2 exp(~k 2 t) with 
h(t) = exp(2 J Q R{t')dt'). The self-consistence condition, 
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(26) 



with the definition tp(k) = (1/VV) J (p(x)d d x, leads to 
the deterministic differential equation for h(t), 



1 ; 

2 



gA 
V 



(27) 



The global magnetization m(t) = <p(0,t)/VV, is just 
given by m{t) — m(0)/ y/h(t). Therefore, m(t) is deter- 
ministic (apart from the randomness of m(0)), and never 
changes sign, which yields 9 = 0. 

Equation ( p7j ) can be solved using Laplace transform, 
but we do not need to know h(t) here. The two times 
correlator is 

A e -k 2 (t+t') 

c(k, t, o = (0(k, tm-K t)) = ; (28) 

The Gaussian process = tp/h(t), has the same corre- 
lator as the diffusion equation. More precisely, <p obeys 
the diffusion equation. Hence the rest of the demonstra- 
tion is the same as above. The persistence exponents 9, 
9q = and the scaling function / are the same as for the 
diffusion equation. 

In these two soluble models, the scaling law of Eq. ( |l6| ) 
is valid for any t and any I, and not only asymptotically 
in the large I large t limit as will be the case in general. 
Remark also that the 9q = result recovers two differ- 
ent behaviors of the global magnetization. For the dif- 
fusion equation, the magnetization is exactly conserved, 
whereas for the large n model it relaxes detcrministically 
to zero. 



we just have to write the evolution equation for the two 
point correlation [fj3| , for t > t', 

2— (r, t, t') = C(r + 1, t, t') + C(r - 1, t, t') - 2C{r, t, t'), 

(29) 

with C(r,t,t') — (S r (t')So(t)}. Summing over r, we get, 

(30) 



2 d(M(t)M(f)) 



at 

hence, (M(t)M(t')) = (M 2 (mm(t, <')))• Then, in the 
scaling regime (M 2 (t)) (x i 1 / 2 , leading to 



Y\ 1/4 

a g (t,t')= [-J , for<><' 



(31) 



The normalized correlator of the global magnetization is 
equal to exp[(it' — u)/4] in the variable u = hit. This 
proves that the Gaussian process M(u) is stationary and 
Markovian and that p(u) ~ exp(— |w|/4). In the t vari- 
able we get $o = 1/4 = (A — d/2)/z, since A = 1 and 
z = 2. Remark that while M(t) is Markovian, S(t) at a 
given point is not, neither is it Gaussian, and the com- 
putation of 9 was a real tour de force [g4| which cannot 
be extended to other systems. 

To check the scaling assumption of Eq. (|l6|) with the 
exact value of 9q, we have simulated the Glauber Ising 
model on a 200000 spins chain with block size 1, 21, 41, 
61, 91. Ten samples were averaged to obtain the final 
data. The data collapse with 6q = 1/A is very good and 
the scaling function has the expected behavior: a power 
law divergence with exponent #o at small argument and 
an algebraic decay with exponent 9 at large argument. 



E. Results for global persistence 

Thanks to the scaling form of Eq. (|l(]), it is possible to 
use block scaling to compute 9q numerically. One eval- 
uates Pi(t) for different I, which can be done on a single 
run, and then adjusts #o to obtain the best data col- 
lapse. We present here some numerical results for three 
different models, illustrating the three possible cases: 9qZ 
equal to, bigger than or smaller than A — d/2. We also 
give a direct derivation of the exact result = 1/4 for 
the one dimensional Glauber dynamics (Majumdar et al 
used an interface representation of the dynamics) . Fi- 
nally we show a surprising relation between 9q for the 
one-dimensional XY model with power law initial spa- 
tial correlations and 9 for the diffusion equation. 

One- dimensional Glauber model - ID coarsening is quite 
special, since the critical temperature is zero. For the 
Glauber Ising model, which is exactly soluble, the exact 
computation of 9 was really difficult, while #o is trivial, 
since the global magnetization M(t) is Gaussian at any 
time, and Markovian in the scaling limit. To show it, 
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FIG. 3. Scaling of the block persistence probability pi(t) 
for a one-dimensional spin chain (200000 spins, 10 samples), 
with block size / = 1,21,41,61,91 (from bottom to top in 
the right part of the insert). I = 1 is omitted in the scaling, 
and the data collapse improves as the block size increases. As 
expected p t (t) ~ t~ 9 ° for t < I 2 and pi(t) ~ t~ e for t > I 2 
(with 0o = 1/4 and 9 = 3/8). 
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One- dimensional model A - In one dimension, determin- 
istic and stochastic models are known to lead to different 
growth laws and correlations p5fl . For instance, in the 
one dimensional noiseless model A, domain walls have 
a weak attractive interaction decreasing exponentially 
with the distance, and L(t) oc hit, whereas in the dis- 
crete stochastic Glauber models, walls do not interact 
but perform simple random walks and annihilate when 
they collide, leading to L(t) oc \ft. 

Model A is interesting, because it can be mapped 
on a simple deterministic model of charge aggregation 
In this model, domains of 



and "-" 

phase evolve the following way. At each step, the small- 
est domain I (length l(Io)) is changed sign and merged 
with its two neighbors I\ and I2, to give a domain of 
length + + l(l2) ■ To compute the domain size 
distribution, the sign of the domains can be forgotten 
and one can easily show that no correlations develop in 
the system. The mean-field equations are exact and can 
be solved for the scaling function of the size distribution. 
In this model, the time variable is the minimum length 

Bray et al. have shown that the local persistence expo- 
nent 9 [^8| and the autocorrelation A [^9| have a geomet- 
rical interpretation in this model. For instance, defin- 
ing for each domain the fraction of persistent spins d(I) , 
the new interval obtained in one step of the aggregation 
model has d(I) — d(Ii) + d(l2), and the total fraction of 
persistent spin can be computed in mean field since there 
are no correlations for d as well. The exact results are 
9 = 0.17504588... and A = 0.6006165.... These exponents 
are solutions of implicit nonlinear equations. 

As far as the global persistence exponent is concerned, 
we have to consider explicitly "+" and "-" domains, and 
now it is quite clear that correlations develop through 
the aggregation procedure. A naive but instructive argu- 
ment neglecting correlations leads to the result 9o = 1/2 
which is in contradiction with the bound 9q < 9. 

We consider a discrete lattice leading to integer values 
of /(/). To increase the minimum domain length in the 
system from Iq to Iq + 1, one has to remove n(Jo) do- 
mains. In each coalescence event, Iq spins change sign 
leading to 8M = ±2/o, depending on the sign of the do- 
main. Now, if we assume that for small M, there is no 
correlation between the sign and the length of domains 
the total SM is a Gaussian zero-mean variable with vari- 
ance oc yj u(Io)21q when u(Iq) 3> 1. It is easily shown 
from the mean-field equations that n(lo) oc 1/Iq in the 
scaling (large Iq regime). Therefore, the magnetization 
increment is a Gaussian variable with constant variance, 
and M(Iq) performs a simple random walk, leading to 
9a = 1/2. This naive value of 9o is clearly wrong, being 
larger than 9. Therefore, to obtain 9g one should treat 
the correlations, which does not seem very simple. 

Numerical results for the domain aggregation model 
are presented in fig. [| for a L = 10 6 chain (30 sam- 
ples). The best scaling was obtained for 9q = 0.165. The 
profile of the scaling function in the cross-over region is 



quite different from the stochastic Glauber model. The 
Markovian scaling law would lead to 9a = A — d/2 0.1 
(z = 1 since the time variable is the dynamical length 
scale). Hence, for this model we have 9q > A — d/2. 
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iog,„(i /D 

FIG. 4. Scaling of the block persistence probability ob- 
tained from simulations of the domain aggregation model 
equivalent to the deterministic one-dimensional model A, on 
a chain with 10 6 sites, l = 71, 81, 101, 111, 131, 161, 191, 211. 
Excellent scaling is found for 80 — 0.165. 



Two-dimensional Glauber Ising model - d = 1 is quite 
special because T — is also the critical temperature, 
and persistence can only be defined at T = 0. Now we 
move to the d = 2 Glauber Ising model, for which block 
scaling will lead to a definition of 9 at finite temperature. 
It is also interesting to determine 9q which could only be 
roughly evaluated by the direct method |^2| despite much 
numerical effort. Comparatively, block scaling is a very 
easy and reliable method. We performed simulations on 
a 2000 2 lattice with blocks of linear size 1,5,9,15,19,25, 
and 31. 20 samples were averaged to obtain the final 
data presented in fig. |5|. We find excellent scaling, with 
9o = 0.09. The uncertainty in the data collapse is roughly 
of 1% on 9q. This value of 9$ is compatible with the range 
0.06 - 0.11 found by Cornell and Sire @. The Marko- 
vian value of z9q would be 11/8 = 1.375 (A = 5/4), and 
for this model we have 9qz < d — A/2. 

Two-dimensional Ginzburg- Landau equation - We can 
also simulate the time dependent Ginzburg-Landau equa- 
tion, corresponding to the continuous model A, 



d t tp = V 2 ip + atp{l - tp 2 ). 



(32) 



Starting from an uncorrelated Gaussian initial condi- 
tion, one can solve the equation using a finite differ- 
ences scheme and compute pi (t) for different block sizes. 
Using block scaling, we can determine 9q and 9. For 
both exponents, we find a value somewhat smaller than 
for the Glauber Ising model: 9q = 0.06 < 0.09 and 
9 = 0.20 < 0.22 (the value of 9 has been also computed 
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by Cornell |30|). For 9, the large time decay of pi(t) 
shows significant curvature and the effective exponent 
seems to increase with time. The scaling function can- 
not be superposed with the scaling function of the Ising 
model. Moreover, the fact that both models have differ- 
ent #o shows that the two-times correlations of the global 
magnetization (which solely determine 9o) are different 
in the scaling regime. This suggests that model A could 
be in a different universality class from the Ising model. 
This was also suggested by Rutenberg |3lJ] in a recent 
paper, as he argued that model- dependent anisotropy in 
the correlation function (due for instance to the lattice) 
does not vanish in the scaling regime. 




0.0 0.5 1.0 

log(t/l 2 ) 

FIG. 5. Block persistence at T — obtained from simula- 
tion of the nonconserved Ising model on a 2000 2 lattice, for 
I — 1,5,9,15,19,25, and 31 (from bottom to top in the in- 
sert), piit) decays as t~ 6 ° at early time and as t~ at large 
time. Excellent scaling is then obtained taking Oo = 0.09. 
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FIG. 6. Block persistence at T = obtained from simu- 
lation of the time dependent Ginzburg-Landau equation at 
zero temperature on a 700 2 grid, for I =9, 11, 13, 15, 17 and 
19 (from bottom to top in the insert). The scaling presented 
is obtained taking 6o = 0.06. The large time decay of pi(i) 
corresponds to 9 = 0.2, with significant curvature though. 



One- dimensional XY model - The one-dimensional non- 
conserved XY model is exactly soluble J32|,[l| , and is quite 
special, as for short ranged initial correlation, the struc- 
ture factor does not exhibit conventional scaling, and the 
growth exponent z = 4 in contrast with general results 
for nonconserved vector spin systems. 

The order parameter is a unitary two-dimensional vec- 
tor ip, defined by its phase a(x,t). The equation of mo- 
tion is simply a diffusion equation 



d t a 



(33) 



with a Gaussian initial condition, 



P({a fe (0)}) ex exp ( - £ y«*(0)a-k 

\ k 



(0)) 



(34) 



The correlation function C(r, £1,^2) = (cos(a(r, fi) — 
a(r, £2)) depends on the initial condition. For a finite 
correlation length in the initial condition (3k = (£/2)k 2 , 
C(r,t,t') does not have the usual scaling form Q|. As a 
consequence, the normalized correlator of the magnetiza- 
tion is not a function of t/t' . 

Now, if the initial correlator has a power law decay 
C(r, 0) ~ r -7 / 71 ', the general scaling form is recovered, 
with 



C(r,t,t') = f 



4tit 2 

Vt+TJ \(t + ?y 



(35) 



with f(x) cx x" 1 ^ UA. For 7 > 7r, the spatial correlations 
of the order parameter decay sufficiently fast for the cen- 
tral limit theorem to be valid (see e.g. |33|). Therefore, 
the magnetization is still Gaussian, and the global persis- 
tence probability, now defined as the probability that a 
component of the magnetization has never changed sign, 
is determined by the the normalized correlator 



a g (t,t') = 



(t + t'f 



(36) 



which continuously depends on 7. This is the local corre- 
lator of the diffusion equation in dimension d — (7— n) /it. 
Since the diffusing field is Gaussian, we conclude that 9q 
for the one dimensional XY model with C(r,t = 0) ~ 
r -ih ig equal to 9 for the diffusion equation in dimen- 
sion (7 — tt)/it, a quite surprising result. 

Before moving to finite temperature, let us mention 
that in a recent paper |34j| , Hinrichsen et al. have used 
our method to study persistence for a directed perco- 
lation model. Their data suggest that 9 = 9 . If this 
surprising result is confirmed, it would be interesting to 
understand the exact reason for it. 
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BLOCK PERSISTENCE AT FINITE 
TEMPERATURE 



Although block persistence is useful even at zero tem- 
perature, our main concern remains finite temperature, 
which we discuss now. The main idea is that because 
the correlation length is finite at finite time, the relative 
fluctuations of the block variables vanish as l~ d l 2 when 
the size of the blocks is increased, and therefore the large 
block limit corresponds to no fluctuations at all, i.e. zero 
temperature. In fact this picture is just a naive justifi- 
cation of the renormalization group flow for coarsening. 
One has to single out T = T c where the relative fluctua- 
tions diverge, as the equilibrium magnetization vanishes, 
and for which a different scaling arises. 



A. T < T c 

Let us first consider < T < T c . The difficulty in 
defining a persistence exponent comes from the fact that 
a spin may flip due to thermal fluctuations, leading to an 
exponential decay p(t) ~ exp(— t/r). Indeed, at T = 0, a 
spin flips only when it is crossed by an interface between 
a + and a — domain, whereas at finite temperature, the 
dominating process at late time, when the domains are 
large, is the flip of a spin within a domain due to thermal 
fluctuations. Therefore, at low temperature, it is natural 
from classical kinetics intuition to expect an Arrhcnius 
law r ~ exp(— AS/T), where AS is the energy barrier 
to flip a spin (or a block) within an ordered domain. As 
T — > 0, r diverges and p crosses over to a power law. 

Arrhenius laws are common enough in physics and 
chemistry, and arise each time a fluctuating process has 
to cross a finite barrier. It is useful, though, to work out 
the random process viewpoint, to clearly understand how 
t should behave with I. 

Let us consider a block of linear size I, and spin block 
variables When L(t) is large enough, the system 
can be considered locally at equilibrium inside a domain, 
and, since there are no long-range correlations, (tpi(t)) « 
l d {if) eq and (Aipt) 2 = (tf(t)) - ( w (t)) a « l d {Av) 2 . 
Therefore the relat ive fluctuation of tpi has the scaling 
A W /( W ) cx y/T/F. 

Thus pi(t) is essentially the probability that the sta- 
tionary random process tpi(t) with mean value of order 
l d and fluctuations of the same order has never crossed 
zero. In other words, it is the survival probability of a 
stationary walker X(t) = (jpi(t)— (<fi))/(<pi(t)), with zero 
mean and a mean square fluctuation (X 2 ) = aT/l d , and 
an absorbing boundary at x = 1. 

To simplify, let us assume that X(t) is Gaussian and 
Markovian. Then one can write a simple Langevin equa- 
tion, 



with a Gaussian white noise ij(t) with (rj(t)r](lf)) — 
2aT/l d S(t — if). In exponential time u = e 27 ', the new 
random variable Y(t) = 2jy/uX performs a simple ran- 
dom walk, 



(38) 



where = rj/y/u is a new Gaussian white noise, with 
(£(«)£(«')) = AjaT/l d S(u - u'). Hence, Pi (u) is the sur- 
vival probability of a simple ld random walker with dif- 
fusion coefficient D = 2aT/l d , starting from x = with 
a moving absorbing wall at x(u) — \fu. The survival 
probability is just, 



S(u) = j P(x, u)dx, 



(39) 



where P(x 7 u) is the presence probability of the walker. 
P(x, u) is the solution of the diffusion equation with an 
absorbing boundary condition at x = y/u. When the wall 
motion is much faster than the diffusion of the walker, 
i.e. D <C 1, which corresponds to small T or large I 
(small fluctuations), P(x,u) can be well approximated 
by a Gaussian distribution with a time dependent weight 
S(u) f§, 



P(x,u) 



S{u) 



VIttDi 



(40) 



where S(u) is determined by equating the mass loss 
rate with the flux of mass through the moving wall. 
At large u, S(u) decays with a power law u~P and 
(3 = (4 7 rL>)- 1 / 2 exp(-l/4L>). Since Pl (t) = S(e 2 ^), we 
recover the heuristic Arrhenius law with, 



r = 1/(27/3) - 



'27raT 
~i 2 l d 



exp[l d /(8aT)} 



(41) 



The constant a is a slowly varying function of the equi- 
librium correlation length but does not depend on I for 
large /. 

The important point is that the effective temperature 
entering the Arrhenius law of the spin blocks is cut by 
a factor l d and that r diverges very quickly when I is 
increased, leading to a fast cross-over to the T = be- 
havior. Admittedly, the actual stochastic process (fi(t) 
is certainly non-Markovian. However, for I much big- 
ger than the equilibrium correlation length, it is nearly 
Gaussian from the central limit theorem. Moreover, its 
correlator C(t) = (<pi(t)(pi(0)) — (<pi(0)} 2 can be bounded 
by two Markovian exponential correlators (because there 
is no long range correlation in time at equilibrium), and 
thus the Arrhenius law still holds with proper constants 
inserted (although the power law in the prefactor may be 
modified), from the discussion of sec. O. 

For t <C r, pi (t) is expected to behave in a similar way 
as for T = 0, and we expect 



X(t) = -jX(t) + V (t) 



(37) 



Pi(t) ~ r z8 °f(t/l*)exp[-t/T(l,T)], 



(42) 
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with two different cross-over times. However, in the scal- 
ing limit / — > oo, t diverges much faster than l z , en- 
tailing that the exponential part does not scale. Hence, 
the scaling form of pi(t) should be Eq. ([l6l). Moreover, 
from the universality of the domain wall dynamics for 
T < T c , the scaling function g should be the same as at 
zero temperature, up to an overall temperature depen- 
dent multiplicative factor. As for the scaling function /, 
we have to take into account a temperature dependent 
multiplicative constant in L(t) (see below). 

To illustrate these ideas, we have performed simula- 
tions of the two-dimensional Glauber Ising model at finite 
temperature on a 1000 2 lattice. Figure ^ presents results 
at T = 2T c /3 for blocks of size I = 1,3,5,7,9, 11,13. The 
exponential decay is clearly visible for I — 1 and 1 = 3. 
However, for larger blocks, r is bigger than the simulation 
time, and pi(t) has the T = behavior, with a power law 
decay with exponent 9 fully compatible with the T = 
value (9 = 0.22), for t 3> l 2 , and a power law decay with 
exponent 9q, for t < I 2 , just as expected. Figure || shows 
the scaling function at T = T c /2 (where the approach 
to scaling is faster) obtained with the zero temperature 
value 9o = 0.09, for I = 7,9, 11, 13. The data collapse is 
really excellent. 
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FIG. 8. Pi(t) expressed in scaling form for T = T c /2, and 
block sizes I = 7,9, 11, 13, using the same value for 28o = 0.18 
as in the T — case. Note the similarity with the T = 
scaling function of fig. 1. 
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FIG. 9. Universality of the scaling function for block per- 
sistence. We show the superposition of scaling data corre- 
sponding to T = and T = T c /2. For T = 0, we (double-log) 
plot / = l 20 °pi(t) versus x = t/l 2 , whereas for T = T c /2 we 
plot / = ail 2e °pi(t) versus x = a^tjl 2 , with a± = 1.07 and 
<J2 = 1.26. The excellent superposition of the two scaling 
functions assesses the universality. 



FIG. 7. Pi(t) for the two dimensional Ising model at 
T = 2T c /3, and block sizes I = 1,3,5,7,9,11,13. The ex- 
ponential decay of pi(t) ~ e~ l ^ T is clearly visible for I = 1 
and I = 3, however for I = 5, the exponential regime is al- 
ready repelled at times longer than the simulation time, in 
agreement with the expected fast divergence r ~ exp(aZ 2 /T). 
For I > 5, only the power law zero-temperature like regime is 
to be seen. 



The temperature universality of the scaling function is 
illustrated in fig. ||. We plot the quantity / = l 20 °pi{t) 
versus x = t/l 2 for a set of zero temperature data, and 
/ = a\l 26 °pi{t) versus x = a 2 t/l 2 for T c /2 data, for blocks 
of size 7,9,11,13. The constants a\ and a 2 are the same 
for all sizes, and are adjusted to superpose the two sets 
of data, ai arises from the temperature dependence of 
the prefactor of t 1 ! 2 in L(t) (the natural time variable) 
and a± is the overall temperature dependent multiplica- 
tive discussed above. The superposition obtained is re- 
ally excellent and assesses the expected universality of 
the scaling of block persistence, in a very similar way as 
what is known for the equal-time two-point spin correla- 
tion function 0] . 
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Thus, block scaling leads to a definition of 9 at finite 
temperature as the exponent of the algebraic decay of 
the scaling function f(x). For the two-dimensional Ising 
model, the temperature independence of 9 obtained with 
this method confirms the results obtained with Derrida's 
definition, but the universality is stronger, since the whole 
block persistence scaling function is universal. This uni- 
versality arises from general arguments and should be 
observed for generic systems. 

Note that universality would rather be expressed in 
terms of L(t) than in terms of t. This is especially 
relevant for the three dimensional Glauber Ising model 
with nearest neighbors interactions on the cubic lat- 
tice. Numerical simulations Hfu} at T = lead to 
9 sa 0.17, whereas at finite temperature our method leads 
to 9t>o ~ 0-26 in agreement with results obtained by 
Stauffer Q using Derrida's definition. For this prob- 
lem, seemingly due to lattice effects, L(t) does not grow 
as t 1 / 2 at zero temperature, but as t 33 [p6[ . At finite 
temperature, lattice effects are overcome, and one recov- 
ers the usual growth law. Now if the block scaling func- 
tion g is universal, one should have the same value of 9z 
at any temperature. From numerical results we obtain 
9z = 3.0 x 0.17 = 0.51 at T = and 9z = 0.26 x 2 = 0.52 
at finite T, which actually confirms this universality. 
Note that these values of 9z are in good agreement with 
an approximate continuous theory p4). 

Block persistence is also very useful to study persis- 
tence for the q-state Potts model, as zero temperature 
dynamics show blocking effects at zero temperature on 
the square lattice with nearest neighbor interactions |57| . 
Working at finite temperature is a more satisfactory way 
of overcoming blocking effects than changing the lattice 
type or including next nearest neighbors interactions. 
Derrida used his comparison method to study the q = 7 
Potts model. His data seemed to suggest a temperature 
dependence of 9 ||. 

On the basis of the present work, we would rather ex- 
pect 9 to be independent of T, at least with our definition. 
To address this question, we have performed simulations 
of the q = 7 Potts model at T = T c /3 and T = 2T c /3, on 
a 1000 2 lattice. We have computed pi(t) for I =1, 3, 9, 
11, 15, 19 and 25, where the block variables are defined 
through a majority rule. For both temperature, accept- 
able scaling is obtained for 9q — 0.3, but the scaling is 
not as good as for the Ising model, and would surely 
be improved by using larger block and simulating longer 
times. Moreover, nonscaling transients extend over quite 
a long period of time for T = 2T c /3. The extraction 
of 9 from the decay of p 9 (t) leads to 0(2T c /3) w 0.485 
and 9(T c /3) = 0.47. The discrepancy is not really signif- 
icant compared to numerical uncertainties, and is much 
smaller anyway than for Derrida's data |9j, who found 
0(2T c /3) « 0.55 and 9(T c /3) w 0.4. Hence, 9 does not 
seem to depend on temperature. This is confirmed by 
the comparison of the two scaling functions, which can 
once again be superposed through a global rescaling. 
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FIG. 10. Scaling of block persistence for the 7-state Potts 
model at finite temperature T — T c /3 and T — 2T c /3, from 
simulations on a 1000 2 lattice for blocks of size I =9, 11, 15, 
19 and 25 (13 samples). Data are noisier than for the Ising 
model, but acceptable scaling is obtained with 6a = 0.3. 



Finally, the value of 9 compares well with zero temper- 
ature data obtained by Derrida et al. These authors 
simulated the next nearest neighbor interactions Potts 
model to avoid blocking effects. Data for p(t) showed 
significant curvature, due to the fact that the effective 
z exponent increases with time, and better results were 
obtained for the exponent ip defined as p(t) ~ L(t)^ v ~ 
E(t) v , where E is the energy difference with the fun- 
damental state. These authors found ip = 1.01 for the 
q = 7 Potts model. Assuming tp = 9z with the asymp- 
totic value z = 2, their data lead to 9 « 0.5, in acceptable 
agreement with our results at finite temperature. 



B. T = T c 

The naive kinetic argument giving the scaling of the 
relaxation rate of pi(t) is bound to break down at the 
critical temperature for several reasons. We know from 
explicit renormalization group analysis p8| ] that T c is a 
fixed point for the dynamics. Hence in contrast with 
T < T c , the thermal decay of the persistence probabil- 
ity must scale. Since the equilibrium magnetization is 
zero, one can no longer make a distinction between a 
slow flip mode due to interface motion and a fast flip 
mode due to thermal equilibrium fluctuations within do- 
mains. There are no domain walls in the system, and the 
relevant length scale is the time dependent correlation 
length ~ t 1 ^" . During the dynamics, patches of cor- 
related spins of length appear in the system. These 
large patches have a large life time due to critical slowing 
down. If we consider a block of spins inside one of these 
patches, the typical time required to flip the spin of the 
block should be roughly speaking of the order of the time 
required to relax a fluctuation of wave vector k 2n/l in 
the critical equilibrium state. Hence, one should have 
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Pi{t) ~ exp(— u> c (2tt /l)t) where u> c (k) is the character- 
istic critical relaxation frequency and scales as k Zc [ j38| . 
This leads to an exponential decay pi(t) ~ exp(— at/l Zc ). 
Therefore, as predicted by the renormalization group ar- 
gument, the exponential decay of the persistence proba- 
bility scales and we must have 



Pl (t)~r^ g m/i) = i-^f(t/n 



(43) 



in the scaling limit I 
f(x) ~ x~ e " when x - 



oo with a fixed ratio £/Z, where 



and f(x) 



when x — > oo. 



The exponent O c is the global persistence exponent at T c 
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FIG. 11. Scaling of the block persistence probability at 
T c from simulations of the two-dimensional Ising model on 
a 1000 2 lattice (15 samples) and blocks of size I =5, 9, 13, 17, 
21 and 27. We have taken z — z c = 2.17 and adjusted 8 C . 
The best collapse is obtained for 6 C — 0.28. 

This scaling theory was checked for the two- 
dimensional Glauber Ising model from simulations on a 
1000 2 lattice and blocks of size 1,5,9,13,17,21 and 27. 
Fifteen samples were averaged to obtain the data shown 
in fig. [ll]. Excellent scaling is found with z c = 2.17 
and 8 C = 0.28, which is also in agreement with the value 
of 9 C obtained by fitting the small x power law. The 
scaling function has an exponentially fast decay at large 
argument, as expected. The value obtained for 8 C yields 
9 c z c w 0.607, some 20% bigger than the value found by 
direct determination of the global persistence probability 
p9]|, and the reason for this discrepancy is unclear. 



VI. CONSERVED MODELS 

The amount of certitudes we have for conserved order 
parameter dynamics is by far much smaller than for the 
nonconserved case. The large n limit is quite pathological 
as it exhibits multiscaling, which is not observed in sim- 
ulations. The only exact result is the celebrated Lifshitz- 
Slyozov- Wagner (LSW) theory for the limit of a 

vanishing concentration of minority phase. In this limit, 



well separated droplets of minority phase are embedded 
in a matrix of majority phase. This spatial structure is 
very different from the labyrinth-like domain structure 
of the equal concentration case. The typical length scale 
L(t) scales as i 1 / 3 (z=3) and one can compute the scaling 
function for equal time correlations. Recently, Lee and 
Rutenberg have shown that A = d for LSW. 

For finite concentrations of the minority phase, and 
especially for the zero-magnetization case, the situation 
was more controversial. While the i 1 / 3 growth law seems 
well established since the numerical work of Husc p2| , 
no conclusive result for A is available. Numerical simu- 
lations of conserved models are difficult because the dy- 
namics are much slower than for the nonconserved mod- 
els, and that corrections to scaling are important even at 
long simulation times (see below). Moreover, the spin- 
exchange Kawasaki dynamics freeze at zero temperature, 
and simulations must be performed at finite temperature, 
and the standard definition of persistence cannot be used. 
This explains why results for conserved persistence can 
be reduced to an analytical computation of 6 in the LSW 
theory pi] ]. In fact, in the absence of numerical simula- 
tions, the question of whether the persistence probability 
has a power law decay or not is open, even if the answer 
is intuitively yes. 

Here, using block persistence, information can be ex- 
tracted from finite temperature data, and we are able to 
study persistence for the Kawasaki d = 2 model, with the 
numerical limitations discussed above. Basically, the dis- 
cussion for nonconserved models can be directly adapted 
to conserved models, but for an important subtleness 
that we now point out. 



A. Block scaling 

In the discussion of the scaling of the block persistence 
probability, some additional care is required compared to 
the nonconserved case. To understand why, let us recall 
the discussion of Majumdar and Huse |2J| for A. 

In the scaling regime, the Fourier transform of the two 
points correlator must have the scaling form C(k, t, t') = 
L d f(kL,L/L'). From the definition of A, we have 
f(x,y) ~ cy~ A at large y. Since the global magneti- 
zation is constant, C(k = 0, i, t') is a time-independent 
constant A (its value depends on the initial condition), 
leading to f(0, y) = Ay~ d for all y, which imposes A — d, 
precisely the result found by Lee and Rutenberg for the 
LSW theory 41 1. However, Majumdar and Huse |2|] 
have shown numerically that the one dimensional con- 
tinuous T = model B, which can be mapped on a de- 
terministic domain aggregation model, has A < d. Why 
does the argument above fail for this model ? 

In fact, we have abusively ruled out the possibility that 
f(0,L/L') = 0. In this case, the k = mode does not 
scale and we must write more generally J25| 



C(k, t, t') = L d f(kL, L/L 1 ) + h(kL, L/L'), (44) 
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where the second term is the leading correction to scal- 
ing. If f(0,L/L r ) = 0, this term is negligible for any 
finite k, but not for k = 0, and we have C(k = 0, t, t') = 
A = fi(0,L/L'). In this case, we obtain no information 
on A, which may be either equal to d or nontrivial. Con- 
versely, if /(0, L/L') is finite, then the previous argument 
applies, and A = d. 

Coming back to block persistence, we make the scaling 
assumption of Eq. ([l6]) . If we just import the result of the 
nonconserved case, we find a — z9q — 0, and f(x) — > 1 
when x — > 0, because of the conservation law. However, 
this result does not always hold. 

The coarse-grained correlator Ci(t,t') is given by, 



c l {t,t')= f -^c(k,t,o 

JkKl ( 27r ) 

In the scaling regime, we have 
d d u 



(45) 



Q{t,t') = 



u<(l/L) 



(f{u,L/L') + L~ d fi(u,L/L')) 

(46) 



Now, in the scaling limit for the persistence probability 
of I — ► oo with fixed ratio x ~ L/l and y = L'/l, we get, 



1 jx 



Ci(t,t')=K d u d - 1 f(u,x/y)du 



(47) 



Now, if we consider the rescaled times x and y, the small 
argument behavior is determined by the small x and y 
asymptotics of the normalized correlator, since in this 
limit, corresponding to L and L' much smaller than /, 
the block variables are Gaussian. If the k = mode 
scales, f(u,x/y) — » A(y/x) d when u — > 0, leading to 
ai(x, y) — ► 1. If /(0, x/y) — 0, we obtain, 



ai(x,y) = n(x/y), 



(48) 



with f(u, x/y) = u 3 n{x/y) + o(u J ). 

Hence ai{x, y) is a nontrivial function of x/y, leading to 
the fact that pi(t) scales as cqx~ 6 , where 9' is a nontriv- 
ial exponent. It means that the constant magnetization 
regime exists at early time, but is confined up to a non- 
scaling cross-over time, and therefore does not appear in 
the scaling function. Finally, the correct scaling for pi (t) , 



P i(t)~r ze f[L(t)/i], 



(49) 



with f(x) oc x~ zS when x — ► oo, and f{x) — x~ z0 when 
x — > 0. The exponent 9' is zero if the k = mode scales, 
and is nontrivial otherwise. Note that 9' can be bigger 
than 9 and that it is possible, although highly unlikely, 
that n(x/y) = (x/y) d+3 \ which would lead to 9' = 0. 

To check the validity of this analytical discussion, we 
have simulated two different one-dimensional conserved 
models, illustrating the two possible cases A = d and 
A < d. 



B. ID-Kawasaki dynamics 



The one dimensional spin-exchange dynamics 
(Kawasaki dynamics) is special, in the sense that it does 
not coarsen at any temperature. Indeed, since its critical 
temperature is zero, coarsening does not occur at finite 
temperature, whereas the system freezes at zero tem- 
perature. However, Majumdar et al. |^3|,^5| have shown 
that coarsening occurs in the T — * limit in the rescaled 
time t = t exp(— 4 J/k^T). The obtained dynamics is 
equivalent to a domain diffusion model of Cornell et al 
J44[ . In this model, domains of length L perform random 
walks with a diffusion constant proportional to X/L and 
coalesce. At small finite temperature, this corresponds 
to the fact that a domain of + phase moves through the 
diffusion of an isolated — spin detached with probability 
exp(— AJ/kbT) from a neighboring — domain and reach- 
ing the other neighboring domain after about L 2 steps 
of a random walk Majumdar et al. have argued 

that for this model A = d, which they have checked 



numerically 43 1. 

The local persistence exponent of the one-dimensional 
Kawasaki dynamics can be defined through this domain 
model. We present on fig. |l2| results of simulations of 
the model on a L — 10 6 chain (10 samples). We ob- 
serve a power law decay with p{r) oc r~ e , with 9 — 0.73. 
To our knowledge, this is the first numerical demonstra- 
tion of the existence of a persistence exponent for a con- 
served model, confirming the result obtained in the LSW 
limit by Lee and Rutenberg |4]J| . Note that the persis- 
tence exponent is much bigger than for Glauber dynamics 
(6 = 3/8) (see below). The complexity of the aggregation 
model leaves little hope of obtaining the exact value of 9 
for this model. 




FIG. 12. Numerical results for the fraction of persistent 
spins in the domain model corresponding to the T — > limit 
of the Kawasaki Ising chain (10 6 sites, 12 samples). We find 
a power law decay at large time with 6 = 0.73. 
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Since A = d, we expect the naive scaling of block per- 
sistence with 9' = 0. We have computed pi(t) for I =101, 
131, 161, 191 and 211, and the persistence probability 
is presented in scaling form with 9' = in fig. [H?]. For 
increasing I, the approach to scaling is very slow and al- 
though we have used very large blocks, the data collapse 
is poor. However, it is worse if one tries to set 9' ^ 0, 
and one clearly sees that f(x) tends to a constant when 
x — > 0, confirming the theoretical prediction. 



o 
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simulations of the domain model on a chain of 10 6 sites 
(20 samples), for blocks of size 101, 131, 161, 191 and 
211. The scaling function presented in fig. [l4| is qualita- 
tively very different from the one in fig. [12]. In agreement 
with the general discussion above and Eq. (|4^), we find 
a cross-over between two exponents 9' > and 9. The 
best data collapse is obtained for 9' = 1.3, whereas we 
find 9 = 0.62 (remark however that the large x decay 
shows curvature). It may seem surprising that 9' > 9. 
However, since 9' is not directly related to the the global 
persistence probability, there is a priori no reason why 9' 
should be less than 9. 




FIG. 13. Block persistence probability in scaling form 
(with 8' = 0) for the domain model corresponding to the 
T — > limit of the Kawasaki Ising chain. Simulations were 
carried out for a 10 6 sites chain (12 samples) for blocks of size 
I =101, 131, 161, 191 and 211 (from top to bottom). The 
approach to scaling is very slow when I is increased, but it 
is clearly visible that the scaling function goes to a constant 
(9' — 0) as expected since A = d. 



C. Deterministic domain model 



FIG. 14. Scaling of the block persistence probability for the 
domain aggregation model equivalent to the ID zero temper- 
ature Cahn-Hilliard dynamics. Simulations were carried out 
on a 10 6 sites chain, for blocks of size 101, 131, 161, 191, 211. 
For this model A < d and the scaling function for block per- 
sistence has a power law divergence at small x. Good scaling 
is obtained with 0' = 1.3, and we find 6 = 0.62. 



Another one-dimensional conserved model, is the zero 
temperature Cahn-Hilliard equation (model B). Majum- 
dar and Huse |25|] have shown that the dynamics could 
be mapped on a deterministic domain aggregation model. 
In one step of the dynamics, the shortest domain Iq of 
length Iq is localized and removed, the left (length h) and 
right (length l r ) neighbors are merged. The length Iq is 
dispatched between the right (l rr ) and left (lu) second 
neighbors (which have the same sign as Iq), according to 
hi = hi + loi and lrr = lrr + Ion with Iq t + Iqi = Iq and 
lor ■ Iqi — h '■ l r - These domain dynamics reflect the fact 
that the shortest domain shrinks due to diffusion fluxes 
from Jo to its second neighbors through its first neigh- 
bors. The fluxes are proportional to e~ 2l °, which makes 
the shortest domain shrink much faster than other do- 
mains. The flux to the right (resp. left) is proportional 
to 1/h (resp. l// r ), which leads to the above ratio of lo r 
and Iqi. 

Majumdar and Huse |^f| have found numerically A = 
0.67. Therefore, this model is in the class A < d and 
we should observe a nontrivial 9'. We have performed 



D. Two-dimensional Kawasaki dynamics 

Effective zero-temperature domain models cannot be 
used successfully to avoid the freezing of the two- 
dimensional Kawasaki dynamics, because of the compli- 
cated geometry of the domains. Using block persistence 
simulations can be performed at finite temperature. As 
mentioned above, simulations are difficult because the 
Kawasaki dynamics are very slow, and do not reach the 
pure i 1 / 3 regime. Therefore, it is difficult to observe block 
scaling and to extract the persistence exponent, and wc 
have to be satisfied with qualitative results. Figure [l5| 
presents data obtained for a 1000 2 systems with a simu- 
lation time of 500000 Monte-Carlo steps for blocks of size 
3, 5, 7, 9, 11, and 15, and 2000 steps for I =15, 21, 25, 35, 
45 and 55. The cross-over in the behavior of pi(t) corre- 
sponding to L(t) ~ I is visible for small blocks. At large 
time, we observe a power law decay with a persistence 
exponent 9 ss 0.5. The actual value of 9 is certainly big- 
ger since the effective z exponent increases with time and 
is still far from its asymptotic value z = 3 (1/z » 0.25 
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at the end of the simulation) . Still, we acknowledge that 
these data are not very conclusive. 

It is not surprising to find a power law decay of the per- 
sistence probability (in the block scaling or T — ► 0) limit, 
because of the slow motion of interfaces, as for the non- 
conserved case. What is less intuitive, is that 9 is much 
bigger for the conserved dynamics than for the noncon- 
served dynamics. However one can understand that fast 
dynamics may lead to a small if one realizes that a 
fast moving domain wall will be ineffective in decreasing 
p(t) if it wipes several times regions of spins that have 
already flipped. Once again, we see that 9 reflects very 
subtle effects. 

At early times, when I is increased, we do not seem to 
have a power law regime, but 9' = 0. According to our 
discussion, this would support the fact that A = d = 2 for 
two-dimensional Kawasaki dynamics, but because L(t) 
has strong corrections to scaling, it is not clear that cor- 
relations correctly scale in the time regime observe, and 
one has to be careful. For large blocks and large time 
with L(t) ~ t 1 / 3, <C I, one might observe a power law. 

10° r- 1 
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FIG. 15. Results of simulations of the Kawasaki 
spin-exchange dynamics on a 700 2 lattice, for blocks of size 
I =1, 3, 5, 7, 9, 11 and 13 (one sample, 500000 steps), 15, 21, 
25, 35, 45, 55 (12 samples, 2000 steps). 



VII. CONCLUSION 

In this article we have introduced the notion of block 
persistence as a generalization of global and local persis- 
tence probabilities, and as a way of giving a meaning to 
the persistence exponent 9 at finite temperature. The- 
oretical arguments as well as results of simulations sug- 
gest that the persistence exponents and the whole scaling 
function of block persistence are temperature indepen- 
dent in the whole T < T c phase, which is conceptually 
speaking very satisfactory. We have also shown that per- 
sistence exponents arise for conserved models as well, and 
that block persistence establishes a distinction between 
two classes of models, A = d and A < d. 

Finally, the important question may be: What do we 
learn from persistence ? In fact, the justification for 



studying such crude models as the one used in coars- 
ening is universality, which states that most of the fine 
details of the system are irrelevant for the study of the 
scaling regime. Then the stake is to identify univer- 
sality classes, to understand the parameters that deter- 
mine them, and also to identify universal quantities. The 
theoretical and numerical study of persistence shows us, 
because it probes temporal correlations very sensitively, 
that model universality is not as wide as it may have been 
hoped a priori from equilibrium-based intuition. While 
universality with respect to initial conditions, or interac- 
tions range, seems to hold in most cases, the present work 
suggests that the continuous model A and the Glauber 
Ising model in two dimensions are in different universality 
classes (at least at zero temperature), even if they have 
the same dimensionality, the same conservation law and 
both short-ranged interactions. Hence, the existence of 
the lattice seems to affect correlations even at large time. 
This could be related to the pertaining of anisotropy at 
large time claimed by Rutenberg in a recent work pl| . 

The authors have benefitted from interesting discus- 
sions with S. Majumdar, B. Derrida, S. Cornell. 
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